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Affine  Frames  of  Rational  Wavelets  in  H*2(n+)  * 
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Abstract 

In  this  paper  we  investigate  frame  decompositions  of  H2(n+)  as  a  method  of 
constructing  rational  approximations  to  nonrational  transfer  functions  in  H2(I1  +  )  . 
The  frames  of  interest  are  generated  from  a  single  analyzing  wavelet.  We  consider 
the  case  in  which  the  analyzing  wavelet  is  rational  and  show  that  by  appropriate 
grouping  of  terms  in  a  wavelet  expansion,  H'(n+)  can  be  decomposed  as  an  infinite 
sum  of  rational  transfer  functions  which  are  related  to  one  another  by  dilation  and 
translation.  Criteria  for  selecting  a  finite  number  of  terms  from  such  an  infinite 
expansion  are  developed  using  time-frequency  localization  properties  of  wavelets. 
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Centers  Program:  NSFD  CDR  8803012,  the  Air  Force  Olfice  of  Scientific  Research  under  contract 
AFOSR-88-0204,  and  by  the  AFOSR  University  Research  Initiative  Program  under  grant  AFOSR-90- 
0105. 


1  Introduction 


Discrete  affine  wavelet  transforms  are  a  class  of  transforms  defined  on  function  spaces 
which  provide  decompositions  of  functions  in  terms  of  dilations  and  translations  of  a  sin¬ 
gle  function  called  the  analyzing  wavelet.  In  this  paper  we  study  wavelet  decompositions 
of  the  Hardy  space  H2(n+)  ,  (where  11+  denotes  the  half-plane  Die  ,?  >  0)  via  rational 
analyzing  wavelets.  A  rational  wavelet  in  H2(H+)  can  of  course  be  viewed  as  the  transfer 
function  of  a  finite-dimensional  causal  linear  time  invariant  system.  We  show  that  any 
function  in  H2(n+)  can  be  written  as  an  infinite  sum  of  real- rational  functions  where 
the  terms  in  the  summation  are  related  to  one  another  by  translations  and  dilations. 
We  refer  to  such  a  decomposition  as  a  Wavelet  System  (WS)  decomposition  of  Id2(n+) 
.  In  requiring  rationality  of  the  analyzing  wavelet,  orthogonality  of  the  ‘basis’  functions 
is  sacrificed.  However,  we  show  that  frames  for  II2(n+)  can  be  constructed  from  dilates 
and  translates  of  rational  analyzing  wavelets.  The  frame  property  assures  completeness 
and  computability  of  the  decomposition.  In  problems  involving  dynamical  systems  with 
nonrational  transfer  functions  (i.e.  infinite-dimensional  systems),  one  is  often  required 
to  approximate  the  system  by  a  simpler  (finite-dimensional)  system.  Equivalently,  one 
is  required  to  contruct  rational  approximations  to  nonrational  transfer  functions.  We 
consider  here  the  problem  of  constructing  a  causal  rational  approximation  to  a  causal 
nonrational  transfer  function  by  selecting  a  finite  number  of  terms  from  the  WS  expan¬ 
sion.  Translations  and  dilations  of  the  analyzing  wavelet  generate  a  family  of  functions 
which  are  well  localized  in  the  time-frequency  plane.  We  discuss  criteria  for  selection 
of  terms  in  a  manner  which  exploits  time-frequency  concentrations  of  affine  wavelets 
and  of  the  original  (nonrational)  transfer  function.  Possible  applications  of  the  methods 
described  in  this  paper  are  in  system  identification  and  compensator  design. 

Central  to  the  theme  of  this  paper  is  the  action  of  the  affine  group  AfF(l)  on  transfer 
functions.  It  is  well-known  that  the  affine  group  Aff(l)  acts  on  transfer  functions  in 
several  ways.  These  were  investigated  in  some  depth  in  [2].  Two  such  actions  are  given 
by  (for  a  >  0,  (3  e  IR), 


(a)  M 

and  (b)  g(s) 


g{as  +  (3) 
<*<7(s) 

1  +  apkg(s) 


The  first  action  can  be  imbedded  in  an  action  of  the  group  S 1  ( 2 )  on  transfer  functions 
due  to  Brockett.  Some  related  group  theoretical  aspects  are  investigated  in  the  paper 
by  Grossmann  and  Morlet  [13].  This  paper  is  closest  in  spirit  to  our  work.  The  key 
difference  is  that  Grossman  and  Morlet  focus  on  continuous  wavelet  transforms  and 
nonrational  wavelets. 


2  Discrete  Wavelet  Transforms  and  Frames  for 
Hilbert  Spaces 

Wavelet  transforms  provide  a  means  for  decomposition  of  function  spaces  via  dilates 
and  translates  of  a  single  function  known  as  the  analyzing  wavelet  (also  referred  to  as 
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the  mother  wavelet).  Although  wavelet  transforms  have  been  studied  in  the  context  of  a 
number  of  different  function  spaces,  a  great  deal  of  attention  has  been  devoted  to  wavelet 
transforms  on  L2(1R)  .  The  contents  ol  this  paper  rely  primarily  on  results  derived  for 
L2(IR)  .  Define,  for  a0  >  0,  the  dilation  operator  Da  on  L2(IR)  by  Daf(x)  -■  a1Hf(ax) 
and  the  translation  operator  Th  on  L2(IR)  by  Tbf(x)  =  f(x  —  b).  Discrete  affine  wavelet 
transforms  on  L2(1R)  allow  representation  of  any  /  6  L2(IR)  as 


/=EE 

m  eZ  71  e Z 


C  V 1 
/7l ,  n  7  771,71? 


where  is  the  analyzing  wavelet  and, 


^m,n(x)  -  Tnb0Damxl){x)  =  a'f 1 2  {(df  x  -  nb0 


(1) 


for  fixed  constants  a0  >  0  and  bQ. 

It  has  been  shown  (see  e.g.  [5])  that  it  is  possible  to  find  functions  if  and  constants  a0 
and  b0  for  which  the  sequence  m,neZ  f°nns  an  orthonormal  basis  for  L2(IR)  .  In 

this  case,  (1)  is  the  usual  orthonormal  basis  expansion  of  a  function,  i.e.  cm>n  =  {/, 
However,  the  requirement  of  orthogonality  of  translates  and  dilates  severely  restricts 
allowable  choices  for  the  analyzing  wavelet.  A  great  deal  of  freedom  in  the  choice  of 
analyzing  wavelets  can  be  purchased  at  the  cost  of  orthogonality  if  one  were  to  consider 
frames  instead. 


2.1  Frames 

Frames,  which  were  first  introduced  by  Duffin  and  Schaeffer  in  [7],  are  natural  general¬ 
izations  of  orthonormal  bases  for  Hilbert  spaces. 

Definition  2.1  Given  a  Hilbert  space  'H  and  a  sequence  of  vectors  {/rn}£L  C  7 i, 
is  called  a  frame  if  there  exist  constants  A  >  0  and  B  <  oo  sxlcIi  that 

^l/ir2<£!</A,  H-’SBiL/I2.  (3) 
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for  every  f  6  hi.  A  and  B  are  called  the  frame  bounds. 

Definition  2.2  Given  a  frame  {hn}  in  a  Hilbert  space  'H,  with  frame  bounds  A  and  B, 
we  can  define  the  frame  operator  S  :  hi  — s-  ‘H  as  follo  ws.  For  any  f  £  'H, 

s/ =£</,*»>/*».  (4) 

n 

The  following  theorem  lists  some  properties  of  the  frame  operator.  Proof  of  these  and 
other  related  properties  of  frames  can  be  found  in  [14]  or  [6]. 

Theorem  2.1 

(1)  S  is  a  bounded  linear  operator  with  AI  <  S  <  BI ,  where  I  is  the  identity  operator 
in  hi. 


•) 


(2)  S  is  an  invertible  operator  with  B~]  I  <  S  1  <  A  '/. 

(3)  Since  AI  <  S  <  BI  implies  that  ||7  —  j^g5||  <  1,  .S’-1 
Neumann  series, 


s-1 


2 

A  +  B 


can  be  computed  via  the 


k 


(5) 


(4)  The  sequence  {S~lhn}  is  also  a  frame,  called  the  dual  frame,  with  frame  bounds 
B~l  and  A~1 . 

(5)  Given  any  f  €  TL,  f  can  be  written  in  terms  of  frame  elements  as 

f  =  Y^<  >  kn  =  <  f'  kn  >  (6) 

(6)  If  there  exists  another  sequence  of  coefficients  {«„}  such  that  f  =  'ffanhn) 

£  kl2  =  £  I  <  /,  s-' k  >  P  +  £  I  <  /,  s-1  K  >  0) 


2.2  Frames  of  Affine  Wavelets 

Given  a  function  in  L2(IR)  satisfying  the  admissibility  condition 


duj  <  oc, 


where  is  the  Fourier  transform  of  i/’,  under  mild  additional  hypothesis,  it  is  possible  to 
choose  parameters  a0  and  b0  such  that  (7 />,«0,  60)  generates  an  affine  frame  for  L2(1R)  , 
i.e.  the  sequence  {i/Vn}  as  defined  by  (2)  forms  a  frame  for  L2(1R)  .  Numerically,  a0 
and  b0  can  be  determined  by  application  of  a  theorem  of  Daubechies  [6]  (see  also  [17]). 


2.2.1  Time-Frequency  Localization  in  Affine  Wavelets 

The  key  property  of  affine  wavelet  decompositions  which  has  demonstrated  noted  ben¬ 
efits  in  a  number  of  applications,  is  the  time-frequency  localization  which  arises  due  to 
the  translations  and  dilations.  The  concepts  of  time  and  frequency  concentrations  of 
a  function  can  be  precisely  defined  (c.f.  [17]),  but  roughly  speaking  these  are  intervals 
in  the  time  or  frequency  domains  which  contain  ‘most’  of  the  energy  in  the  signal.  Let 
=  [u>o(ip),  denote  the  frequency  concentration  of  f  and  R(if)  =  [t0(f>), 

(t0  >  0)  denote  the  time  concentration  of  (/',  where  f  is  an  admissible  analyzing  wavelet. 
Then  </>  is  a  function  which  is  concentrated  in  the  time  frequency  plane  on  the  rectangle 
Q  =  O(V0  X  R(ip)  and  each  of  the  wavelets  are  concentrated  on  rectangles 

Qm,n  =  +  nho),%m{tA'l>)  +  nbo)]  X  [«™^o(^).«mWi(^)]. 

(Assuming  symmetry  about  u  =  0,  only  positive  frequencies  are  considered.) 
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3  Wavelet  Decompositions  of  ID (11+) 

3.1  The  Hardy  Space  H2(I1+) 

The  class  of  transfer  functions  of  interest  to  us  is  the  Hardy  space  H2(fl+)  where 
n+denotes,  the  half-plane  9ie  .s  >  0.  For  completeness,  a  few  important  properties 
of  H2(n+)  are  reviewed  here.  For  further  discussion  on  IT’  spaces  see  [15,  8,  9].  In 
particular  [15]  discusses  the  space  kT(n+)  and  the  relationship  between  Hp  spaces  on 
the  unit  disc  and  Hp  spaces  of  11+ 

Definition  3.1  Given  a  function  F  which  is  analytic  in  11+ ,  F  is  said  to  belong  to  the 
class  H2(I1+)  if 

/oo 

\F(x  +  iy)\2  dy  <  oo.  (8) 

-oo 

H2(II+)  is  a  Banach  space  with  norm  (denoted  ||  •  |j/./2 )  defined  by  (8). 

Some  of  the  most  basic  properties  of  H2(n+)  are  captured  by  the  following  theorem 

(c-f.  [15]). 

Theorem  3.1  Given  F  G  H2(Il+);  the  following  are  true: 

(1)  The  nontangential  limit  of  F  exists  at  almost  every  point  on  the  imaginary  axis. 
(3)  The  boundary  value  function  of  F  is  in  L2(IR)  and, 


F(x  +  iy)  =  -  f  F{iu>)- 
Jr 


-du,  x  >  0 


\ w  if  f  I  - /  -2  ,  /  \-)  ' ’  w  v 

7T  Jn  x2  +  (y-u)- 

(3)  The  functions  Fx(y)  —  F[x  +  iy)  converge  in  L 2  norm  to  F(iy)  as  x  —■  0. 

The  elements  of  H2(I1+)  are  transfer  functions  of  causal  input-output  stable  linear  sys¬ 
tems.  More  precisely,  we  have  the  classical  result, 

Theorem  3.2  (Paley-Wiener)  /l  complex-valued  function  F  is  in  Il2(n+)  if  and  only 

if, 

r  oo 

F(s)  —  /  f{l)e~st(il, 

Jo 

for  some  f  G  L2(0,oo)  .  Furthermore  this  representation  is  unique. 

By  the  Paley-Wiener  theorem, 


- —  [  F(iu>)e,wtduj  =  0  for  t  <  0. 
Jr 


Hence  boundary  values  of  functions  in  Il2(n+)  comprise  a  subspace  V  of 
L2(IR)  characterized  by  (9).  We  also  know  by  Theorem  3.1  that  given  the  boundary 
value  function  F(iui),  F  can  be  recovered  on  the  right  half-plane  by  the  Poisson  integral. 
Also  ||P||W2  =  \\F\\L2  where  F  is  the  nontangential  limit  of  F. 

A  Hilbert  space  structure  is  placed  on  H2(ll+)  by  the  inner  product  defined  by, 

ro © 

(F,G)h  2—  /  F(iu>)G(iu>)du. 
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Remark:  In  what  follows,  we  take  the  liberty  of  identifying  H2(FI+)  with  V  .  We  use 
u  to  denote  the  real  frequency  variable  and  write  F{u;)  to  denote  the  boundary  value 
function  F(iu >).  The  complex  frequency  variable  will  be  denoted  by  s.  In  the  remainder 
of  this  paper,  we  adopt  the  following  notation: 

RH2(n+)  :  Real-rational  functions  in  H2(fl+)  . 

H2  (11+)  :  Functions  in  H2(n+)  which  are  Laplace  transforms  of  real-valued  functions  in 
L2(0,  oo)  . 

Note  that  RH2(n+)  C  H2l(n+). 


3.2  Wavelet  Transforms  on  H2(n+) 

By  the  Paley-Wiener  theorem  RH2(H+)  represents  transfer  functions  of  causal  finite¬ 
dimensional  linear  systems  with  square  integrable  weighting  patterns.  Of  particular- 
interest  to  us  is  the  case  where  the  analyzing  wavelet  T  belongs  to  the  class  RH2(LI+). 
For  any  F  E  H2(II+)  define  the  restriction  of  F  to  vertical  lines  in  the  right  half-plane 

by, 

Fx[y)  =  F{x  +  iy),  X  >  0. 

Also  define  the  Fourier  transform  along  vertical  lines  in  II+by, 

Friu)  =  7T  I  F(x  +  iy)e‘uydy 

dip 

for  F  E  H2(fl+). 

Now  let  T  e  RH2(II+)  C  H2(II+)  be  our  candidate  for  an  analyzing  wavelet,  and 
suppose  that  for  any  r>  0  $  satisfies  the  admissibility  condition, 


CV 


,M“) 


'r  M 


-du  <  oo, 


x  >  0. 


(10) 


Under  these  assumptions, a  family  of  continuous  wavelet  transforms  can  be  defined  on 
H2(n+)  .  Let, 


TUU(S)  =  | a|1/2  T(os  -  ib),  a,  b  E  IR,  a  >  0. 


Then  for  any  F  E  H2(I1+)  define  the  continuous  wavelet  transform  on  the  line  3?e  s  =  x 

by, 

r 

(11) 


WxF(a,  b)  —  /  Fx(y)^:'bKy)dy. 


Inversion  of  this  transform  is  accomplished  by  the  inversion  formula, 


FM  = 


1 


a 


Wx.F(a,  b)  vkU,6)(:f/ ) 


da  db 


JRJlt 


To  define  a  discrete  wavelet  transform  with  respect  to  T,  let 

'km,n(-s)  =  aul2\])(a^s  -  inb0)  «0  >  0.  3fe  s  =  x. 


Assume  T  satisfies  the  admissibility  condition  and  let  a0  >  0  and  b0  be  such  that  the 
family  {^ra,n}m„eZ)  is  a  frame  for  V  which  is  a  closed  subspace  ot  L2(IR)  .  Thus  for 
any  F  €  H2(fT+)  we  have  that 

m  71 


where  S  is  the  frame  operator  associated  with  the  frame  { vhmn }  • 

4  Wavelet  System  Transfer  Functions 

In  this  paper  we  restrict  discussion  to  affine  frames  constructed  on  the  imaginary  axis. 
Hence  we  require  admissibility  of  the  analyzing  wavelet  T ( s )  for  3?e  s  =  0  and  consider 
the  boundary  values  of  Id2  (11+)  . 

Let  ®(w)  denote  the  nontangential  limit  of  an  admissible  analyzing  wavelet  'l(s) 
and  let  (a0,  b0 )  be  such  that  ( ,  «0,  b0 )  generates  an  affine  frame  for  Ii2(n+)  .  Then  any 
F  6  H2(n+)  can  be  represented  as 

F{w)  =  X]  S  (12) 

m  n 


where  as  before  —  nb0).  Assuming  that  F  £  IT2(n+)  is  the 

Laplace  transform  of  a  real-valued  weighting  pattern  in  L2(0,  oo)  ,  and  $  is  a  real-rational 
wavelet,  we  can  ask  the  question;  is  an  arbitrary  truncation  of  the  frame  expansion  (12) 
the  transfer  function  of  a  real  weighting  pattern?  Since  RH2(n+)  arises  as  the  image 
under  the  Laplace  transform  of  certain  real-valued  weighting  patterns,  if  the  analyzing 
wavelet  $  e  RH2(n+)  then  Tmn(-cu)  =  T*nn(u;)  for  n  ~  0  and  ni  e  Z.  However, 
this  symmetry  is  violated  for  n  y  0.  To  circumvent  this  problem, we  can  group  together 
positive  and  negative  translates  of  the  analyzing  wavelet  to  form  what  we  call  a  wavelet 
system  transfer  function. 

Definition  4.1  Given  G  RH2(n+)  C  H2(n+)  a  wavelet  system  transfer  function  is 
defined  as  a  function  of  the  form, 

Gm’n(s)  =  a\Sfmtn(s)  +  m imt_n(s) 

where  a  denotes  the  complex  conjugate  of  a. 

Proposition  4.1  Let  F(s)  be  a  real-rational  function  in  H2(n+)  .  Then,  the  wavelet 
system  transfer  function  defiiied  by 

Gm’n(s)  =  aFm<n(s)  +  aFm_n(s), 

is  also  a  real- rational  function  in  H2(n+)  . 

Proof:  Let  F(uj)  and  denote  the  nontangential  limits  ol  F{s)  and  Gm’n(s) 

respectively.  Since  H2(n+)  is  a  closed  linear  space,  it  is  clear  that  Gm'n(s)  £  H2(n+).  It 
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is  also  clear  that  Gm'n(s ) 
to  show  that 


is  rational.  To  show  that  G'm>n(.s) 

/  \  -prill, n  .  > 

Gm’n{-uj)  =  G  (t o). 


is  real-rational  we  only  need 


Now, 


Gm,n(—u>)  = 


0-Fmr i(  T  Oil  m  _ n[  to) 

aa^“F{— «”lto  -  r?i>0)  +  aao"^2F(— «”lco  +  nb0) 
aa™/2F(-(a™co  +  nb0))  +  aa™12 F(-{a™ to  -  nb0 )) 
aa^2  F  (a™  u  +  nbQ)  +  cx(Iq/2  F(a"lu  -  nb0) 
o^Fm  ^ ?1(to)  T  aFmn(iu) 


Grouping  together  translates  of  the  analyzing  wavelet  to  form  WS  transfer  functions 
as  in  Definition  4.1,  requires  that  the  coefficients  that  appear  in  front  of  each  of  the 
two  individual  wavelet  terms  be  complex  conjugates  of  one  another.  For  WS  transfer 
functions  to  be  useful,  we  need  to  verify  that  the  coefficients  (A,  S~l  tym<n)  satisfy  this 
condition  for  a  interesting  class  of  functions  in  H-(H+)  .  In  verifying  this  condition,  we 
make  use  of  the  following  proposition. 

Proposition  4.2  Let  f  and  g  be  boundary  values  of  functions  in  H2(fl+)  .  Also  let  f 
and  g  be  such  that  /(— to)  =  /(to),  and  g(—  to)  =  g( to),  then, 


^\fi9m,n)  (,J  i  9m,— n)  ■ 


Proof:  (See  Appendix  A) 

Using  Proposition  4.2,  we  can  show  that  for  functions  in  IT2(n+)  ,  which  are  the  Laplace 
transforms  of  real-valued  functions  in  L2(0,  oo)  ,  the  wavelet  coefficients  do  indeed  satisfy 
the  conjugation  condition  required  for  utility  of  the  concept  of  WS  transfer  functions. 

Lemma  4.1  Let  f  and  g  be  boundary  values  of  functions  in  H2(n+)  such  that  /(—to)  = 
/(to),  and  g(— to)  =  g(w).  Let  a0  >  0  and  b0  be  such  that  { a  frame  and  let 
S  be  the  associated  frame  operator.  Then, 


Om,n  (/)  S  9m,n )  (7) ‘5  9m,— n)  G 


Proof:  Since  the  frame  operator  S  is  self-adjoint, 

9m, n )  {*5  1 1  9m,n )  • 


Therefore  by  Proposition  4.2  we  need  only  check  that  if  /(—  to)  =  /(to)  then  S  1  /(—to)  = 
S~lf{ to).  Now, 


(s-7)(-w) 


2 

A+  B 


A  +  B 


S 


/(-to). 


(13) 


i 


Let 


Gm'n(u) 

Gm'°{u) 

Then 


Thus, 


(/ ■>  9m,n)  9m,niW)  T  (./)  SW,— n)  9  m,— ni^)  D1  6  Z  ,  11  £  Z  ,  11  -f-  0 

(/>  Sm.o)  Sm.oM  »>eZ. 


(s/k-w)  =  £ 

m 

-  E 


Gm-°(-w)  +  ^  G'm- 


n=l 

oo 


Gm’°(u>)  +  J3  G'm,n(u;) 


71  =  1 


(5/)(«). 


>4  +  B 


5  / 


(-«)  =  /- 


/!  +  £ 


5  /M 


and  together  with  (13),  we  get  (,S~l —  (S~l which  completes  the  proof. 


The  above  observations  lead  us  to  the  definition  of  what  we  will  call  the  Wavelet 
System  Transform  on  H2(IT+)  .  We  first  define  the  space  containing  the  image  of  the 
WS  transform  operator. 

Definition  4.2  Define  7 Z(J)  as  the  space  of  sequences  of  functions  {Fk}kej  su°h  that 
Ff.  £  RH2(n+)  C  H2(n+)  for  all  hej  and 

II  53  F^ln  <  00 •  (14) 

fcej 

Definition  4.3  Let  'L  £  H2(II+)  be  an  admissible  real-rational  analyzing  wavelet  and 
a0  >  0,  b0  such  that  (ty,a0,b0)  generates  an  affine  frame  for  H2(fl+)  .  The  Wavelet 
System  transform  is  defined  by  the  operator  H $  :  H|l(n+)  — >  7 Z(7L  x  7ZF), 

mF={F^)M z  „6Z+. 

where  for  any  F  E  H2(L1+),  Fm<n  is  given  by, 

Fm’n(  a)  =  +  (F,S_14m>_n)  <t„..n(w)  me  Z,neZ+,  re  ^  0 

F”-°M  =  ire)  m€  Z.  (15) 


The  operator  By  defines  an  invertible  isometry  from  H2(TI+)  to  7Z(7L.  x  2+),  (where  (14) 
defines  the  norm  on  Inversion  of  //,(,  involves  a  simple  summation  of  terms  in 

the  sequence  { Fm'n }. 

Putting  together  Proposition  4.1  and  Lemma  4.1,  we  can  now  state  the  main  decom¬ 
position  theorem  for  H^(n+)  . 
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Theorem  4.1  Let  'X'  £  RH2(ff+)  be  an  admissible  analyzing  wavelet,  and  a0  >  0,  bQ 
such  that  (\P,a0,60)  generates  an  affine  frame  for  II 2 ( T1 + )  .  Then,  any  F  in  H'j^(n+) 
can  be  represented  as, 

CO 

f  =  EEf*'  (16> 

m  n= 0 

where,  each  Fm’n  (e  RH2(II+)y)  is  a  wavelet  system  transfer  function  defined  by, 

pm.n  -  [F,  S~1]illhn)  +  (F,S~^'ilm  n)'imi_n  m  G  Z,  n  =  1,2.... 

Fm'°  =  {F,S~^ PTO,0>  ^m,0  me  Z.  (17) 

Note  that  if  the  analyzing  wavelet  If  is  of  degree  N ,  then  the  degree  of  each  Fm'n  is 
bounded  by  2 N. 


4.1  Example:  Rational  Analyzing  Wavelet  for  H2(I1+) 

As  an  example  of  a  rational  analyzing  wavelet  for  H2(fl+)  ,  consider  the  function 


'l(.s) 


1 

(s+7)2+r 


7  A  >  0. 


It  is  easily  verified  that  \P  (s)  is  in  H2(II+)  and  furthermore  'P  is  an  admissible  analyzing 
wavelet  i.e., 


/  %{x  +  h 

J  R 


y)dy  =  0 


for  x  >  0. 


Taking  the  inverse  Laplace  transform  of  \P ,  the  weighting  pattern  if(t)  corresponding  to 
4/  is  given  by, 


fi(t)  = 


£— ie-7fsin (ft)  for  t.  >  0 
0  for  t  <  0 


The  H2(II+)  norm  of  *P  can  also  be  found  to  be, 


27r  / 1  cos(arctan(£/7)) 

H 2  W  \7  VC  +  7 2 

Figure  1  and  2  show  this  wavelet  'P  evaluated  on  the  imaginary  axis,  and  the  corre¬ 
sponding  weighting  pattern  if  for  different  values  of  f  and  7. 


4.2  Constructing  an  Affine  Frame  for  H'2(II+)  from  'F 

For  the  purpose  of  numerically  determining  values  of  the  dilation  stepsize  a0  and  trans¬ 
lation  stepsize  b0  such  that  (IP,  a0,  b0)  generates  an  affine  frame  for  H2(II+)  we  can  utilize 
Theorem  C.l  and  Corollary  C.l  in  Appendix  C.  Figures  3-4  show  the  results  of  applying 
Theorem  C.l  and  Corollary  C.l  for  the  case  where  a0  —  2.0.  Hence  for  a0  =  2.0  and 
0  <  b0  <  17  the  sequence  {vp„lin}  forms  an  affine  frame  for  H2(n+)  . 
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w  t 

Figure  1:  H2(II+)  wavelet  T  for  7  =  5,£  =  1,  and  weighting  pattern  ij>. 

5  Properties  of  Wavelet  Systems 

5.1  Time-Frequency  Localization 

Perhaps  the  most  useful  property  of  wavelet  system  decomposition  is  the  time-frequency 
localization  which  arises  due  to  the  translations  and  dilations.  As  in  Section  2.2.1, 
let  O('I')  =  [o;0(vl/),c<;1('l)]  denote  the  frequency  concentration  of  T  and  R.('b)  = 
[t0(df),  f1(d')]  (t0  >  0)  denote  the  time  concentration  of  T  (computed  using  the  inverse 
Laplace  transform).  Then  'L  is  a  function  which  is  concentrated  in  the  time  frequency 
plane  on  the  rectangle  Q  =  fl(d')  x  /?( VP )  and  each  of  the  WS  transfer  functions  are 
concentrated  on  rectangles 

Qmn  =  [aomM'S>)  +  nbQ),a-m(u1('lf)  +  nb0)}x[a”lt0(<!!),a'nt1('S>)}  =  fl(d 'm,„)  x /2($m,n). 

Figure  5,  shows  the  distribution  of  the  rectangles  Qmn  in  the  time-frequency  plane. 
Thus  the  WS  expansion  (16)  provides  a  decomposition  of  H2(II+)  in  which  each  term 
represents  a  finite-dimensional  system  which  is  localized  in  both  time  and  frequency. 
Since  the  wavelets  for  H2(II+)  are  constructed  in  the  frequency  domain,  the  time  and 
frequency  axes  are  interchanged  when  compared  to  the  analogous  picture  for  wavelets 
constructed  in  the  time-domain.  A  noteworthy  feature  of  time-frequency  localization 
properties  of  H2(n+)  wavelet  systems  is  that: 

•  Near  t  =  0,  while  time  localization  is  good,  the  frequency  concentration  of  each 
wavelet  system  encompasses  a  large  band  of  frequencies. 

•  For  t  0,  frequency  localization  is  very  good  so  that  it  is  possible  to  ‘zoom’  in  on 
narrow  frequency  bands,  whereas  localization  in  time  is  poor. 
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Figure  2:  H2(n+)  wavelet  $  for  7  =  5,£  =  10,  and  weighting  pattern  0. 


5.2  Poles  and  Zeros 

Given  the  poles  and  zeros  of  the  analyzing  wavelet  'I'  we  would  like  to  examine  the 
manner  in  which  the  poles  and  zeros  of  the  WS  transfer  functions  are  influenced  by  the 
translations  and  dilations.  Assume  first  of  all  that  T  has  no  real  axis  poles  and  that 
is  a  degree  N  real-rational  function  in  H2(fl+)  .  Let  {PkiPk)k=i  '■)e  die  se^  set  P°^es  °f 
'F(s),  and  {zj,z ~}J=1  be  the  set  of  zeros  of  T.  Since  T  is  real-rational,  'L  can  be  written 


*{*)  = 


P(s)  ni(*-*i)(*-**) 


w  Q(s)  n  k(*-pk)(*-pi;Y 

where,  P(s)  and  Q(s)  are  relatively  coprime.  Thus, 


= 


ao  ~  n }(aos  ~  ZJ  -  inb0){a™s  -  -  inbQ) 

ILK**  -  Pk  -  ™bQ)(a™s  -pH-  inb0 ) 

<hT/2  n,-(^  -  ftj("M*))(3  -  Kj(m,n))  _  PnAs) 

U))(S  -  "d™,  n))  Qm,n(s)  ’ 


where. 


/3f.(m,n)  =  aom(zj  +  inb0)  7  k(m,n)  =  aQm(zj  +  inb0)  (19) 

r]k{m,n)  =  a~m(pk  +  inb0)  uk(m,n)  =  a~m(pj;  +  inb0)  (20) 

Note  that  T]k(m,—n)  =  T^(m,n)  and  i/k(m,—n)  =  7/i~(m,n)  and  /3k ;(m,  — n)  =  7 k(m,n) 
and  7 k(m,  —  n)  =  (3k(m,n ).  Therefore  the  poles  of  Gm’n(s)  =  or^miTl(s)  +  cfTm  _n(s)  are 

{Vk(m,7i),7]j:{7n,n),  ^(m,n),]^(?ii,n)}  .  (21) 

Thus  the  effect  of  dilations  and  translations  upon  poles  is, 
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Figure  3:  Estimates  of  frame  bounds,  using  H2(T1+)  analyzing  wavelet  T  of  Example  1, 
with  dilation  stepsize  a  —  2,  as  translation  stepsize  b  is  varied.  Solid  curve  represents 
the  lower  frame  bound  A  and  the  dashed  curve  represents  the  upper  frame  bound  B. 

•  Dilations  move  the  poles  of  Gm'n  radially  away  from  and  towards  zero.  As  the 
dilation  index  m  increases,  the  poles  move  towards  zero  and  as  m  decreases  the 
poles  move  away  from  zero.  However,  the  poles  remain  in  the  closed  left  half-plane 
which  is  crucial  since  otherwise  the  Gm<n  would  not  remain  in  H2(H+)  . 

•  The  complex  translations  simply  translate  the  poles  along  vertical  lines  in  the  left 
half-plane. 

If  we  write  Gm’n(s)  =  Nmjl(s)/ Dmn(s),  then  the  zeros  of  Gm<n(s)  are  the  roots  of 

Nm,n0)  =  aPm<n{s)Dm^_n(s) +  aPm_n{s)Dm<n{s) 

=  a  ~  n))(s  “  7j(™> «))  n))(5  “  n)) 

j  k 

+a  -  Pj(m,  n))(s  —  7j(m,  n))  ]j[(s  -  ?/fc(m,n))(s  -  vk (m,n)) 
j  k 

Without  actually  deriving  expressions  for  the  zeros  of  Gm,n(s)  we  make  the  following 
observations: 

•  The  zeros  of  Gm’n(s )  occur  in  complex  conjugate  pairs. 

•  For  n  ^  0,  if  Gm<n  has  no  real-axis  poles,  then  Nm  n(s)  and  Dm  n  are  coprime,  and 
therefore  Gm'n(s )  is  a  strictly  proper  rational  function  of  order  2 N. 

•  For  n  =  0  pole-zero  cancellation  results  in  Gm,n(s)  of  order  N. 

•  If  Z(n)  is  the  number  of  poles  of  Gm'n  on  the  real  axis  then  the  order  of  Gm'n(s)  is 
2  N-l(n). 


70 


60 


Translation  Stepsize  b 

Figure  4:  Ratio  (B/A)  of  estimated  frame  bounds,  using  H'2(II+)  analyzing  wavelet  'I'  of 
Example  1,  with  dilation  stepsize  a  =  2,  as  translation  stepsize  b  is  varied.  Solid  curve 
represents  B/A,  and  the-  dashed  line  indicates  the  level  where  B/A  =  1. 


5.3  State  Space  Realizations  of  Wavelet  System  Transfer  Func¬ 
tions 


Let  if(t)  be  the  weighting  pattern  (i.e.  inverse  Laplace  transform)  corresponding  to 
the  real  rational  H2(IT+)  -wavelet  $  and  define  the  triplet  {A,B,C)  to  be  a  minimal 
realization  of  In  this  section  we  examine  the  question  of  realizations  of  the  WS 
transfer  function  Gm’n{s)  =  a^m,«(5)  +  a*\ _n(s).  Let  aR  and  <*,  denote  respectively 
the  real  and  imaginary  parts  of  o.  Since  (.4,  B,  C )  is  a  realization  of  VP, 

?/>(£)  —  C  exp  (At)B. 


Taking  the  inverse  Laplace  transform,  the  weighting  pattern  corresponding  to  Gm’n(s ) 


is. 


9m,n(t)  =  aa-m/2f(a-mt)e^o'n^t  +  aa-"*/20(a-™i)e-‘noo  m»°‘ 

=  aQm/2if(a-mt)2  [aRcos(na~mb0t)  -  a,  sin(na-m60<)]  (22) 

=  aaom/2Cexp{a-mAt)Beinao,rXbot  +  aafm/2C  exp{a~m At)Be-inaon'b°t 
=  2 C  exp((a~m  A  +  ina~mb0I)t)B 

+aafml2C  exp((a“"M  -  ia~mnb0I)t)B .  (23) 

Equation  (23)  immediately  provides  a  complex  realization  of  Gm’n(s)  which  is 


A  = 

n 


a~mA  +  ina~mb0I 


0 


0 

a.frn  A  —  inafmb0I 


B„ 


=  [B  B\‘ 


r  - 

W  77171 


i  :m/2aC  a:m/‘aC 


(24) 

(25) 

(26) 
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(Translations)  frequency 


Figure  5:  Time- frequency  concentrations  of  wavelet  systems. 


However,  since  V\n,n(0  is  a  rea^  weighting  pattern,  we  need  a  change  of  basis  to  give  us 
a  real  realization.  The  differential  equations  corresponding  to  (26)  are 

a?!  =  (a“m,4  +  ina~mb0I)x1  +  Bu 

x2  =  -  ina~mb0I)x2  +  Bu  (27) 

and  the  output  map  is  given  by 

y  =  a~m/2aCxx  +  «o,n/2cv  Cx2.  (28) 

Let,  zx  =  (aq  +  x2)  and  z2  —  i{xx  —  x2).  Then  we  get 

zx  =  a~mAzx  +  na~mb0z2  +  2Bu 

z2  =  fl-m/l22-n(i-mt0:1.  (29) 

Under  this  transformation  the  output  map  becomes 

V  —  a0m^2  {aRC{xx  +  x2)  +  iaIC(xl  —  x2)) 


=  a0  m/2  ( aRC zx  +  iajC z2)  . 

(30) 

From  (29)  and  (29)  we  get  the  following  real  realization  of  Gm'n(s), 

[  n~~A  ' 

Amn  ~  [  -na^'bj  a-’“A 

i  n  #  0 

(31) 

B,„„  =  t 2B  Of.  n?0 

(32) 

Cmn  =  [«0  ",2<*rC  -1 T/2«,C 

,  n  /  0 

(33) 

(/!„,„,  Bmi„,  <?„,,„)  = 

(34) 

It  is  interesting  to  note  that,  in  this  form  (31-33),  the  dilations  and  translations  which 
appear  in  the  WS  transfer  function  Gm>n(s)  affect  the  state  space  realization  via  dilations 
and  translations  as  well,  i.e. 
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5.3.1  Parallel  Connections  of  WS  Transfer  Functions 


In  constructing  rational  approximations  (see  Section  7),  we  will  be  considering  parallel 
connections  of  WS  transfer  functions  of  the  form, 


(m,n)&J 


(35) 


where  J  is  a  finite  index  set.  Given  Equations  (31-33)  defining  the  realization 
(Am  n,  Bm  n,  Cm  n)  of  Gm'n(s):  a  state  space  realization  of  G j(s)  is  readily  obtained. 


Aj  - 


A 


771 J  ,7Xj 


0 


■^7711  ,712 


A  2 

1  A  7712  ,71 1 


A  2 


(36) 


Cj  ,n\  1  '  '  J 


(37) 

(38) 


5.3.2  Minimality  of  State  Space  Realizations 


An  obvious  question  that  can  be  asked  regarding  the  realization  (36-38)  of  Gj(s)  is 
whether  such  a  realization  is  minimal  in  the  sense  of  being  both  controllable  and  ob¬ 
servable.  We  have  the  following  result, 

Theorem  5.1  Let  e  H2(II+)  be  an  admissible  rational  analyzing  wavelet,  let  [pj } 
denote  the  poles  o/'h  and  let  ( A,B,C )  be  a  minimal  realization  of  'k  (of  dimension  N ). 
Also  let  Jbe  a  finite,  bounded  index  set.  Then 


(a)  The  realization  (31-33)  of  Gm'n(s )  is  minimal  for  all  (m,n)  E  J  if  and 
only  if,  for  any  j  =  1, . . . ,  N ,  there  do  not  exist  any  nonzero  integers  k,  such 
that. 


Im  pj 


=  k. 


(39) 


(b)  The  realization  (36-38)  of  G j(s)  is  minimal  if, 

( 1 )  {Amn,Bm^n,Cm<n)  as  defined  in  (31-33)  is  a  minimal  realization  of 
Gm'n(s)  for  all  ( m,n )  E  J ,  and 

(2)  For  any  j,l  =  1, . . . ,  N ,  there  do  not  exist  any  nonzero  integers  k,  such 
that. 


1  1  f  Pj\ 

log  «0  °s  V  Pi  J 


=  k. 


(40) 


Proof:  (See  Appendix  B)  ■. 
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6  Computing  the  Wavelet  System  Decomposition 


In  principle,  given  a  frame  the  coefficients  of  a  WS  decomposition  of  any 

f  €  H2(n+)  can  be  computed  by  first  computing  the  dual  frame  {5'~1'I/m  n}  using 
the  Neumann  series  expansion  for  the  inverse  frame  operator  S'-1  and  then  computing 
the  inner  products  of  /  with  elements  of  the  dual  frame.  However  in  practice,  this  may 
prove  to  be  an  extremely  cumbersome  computation,  especially  in  cases  where  the  frame 
is  not  close  to  being  a  tight  frame  (B/A  =  1).  Expansions  with  respect  to  frames  are  in 
general  not  unique  due  to  the  lack  of  linear  independence.  It  can  be  seen  from  Theorem 
2.1  that  among  all  possible  expansion  coefficients  for  a  given  function  /,  the  coefficient 
sequence  c  ={(/,  }  is  optimal  in  the  sense  of  minimum  t2(7L2)  norm.  That  is 

if  { am,n }  is  such  that 


/-EE 


a  T 

^  ?n,7i 


m  n 


then, 

EEK/.5-,^)l2<EEI“...d2- 

ran  rii  n 

For  WS  approximations  using  a  finite  number  of  terms,  the  problem  of  computing  co¬ 
efficients  may  be  posed  as  a  least  squares  approximation  problem.  Among  all  possible 
sequences  of  approximation  coefficients,  the  minimum  t2  norm  coefficient  sequence  can 
be  obtained  by  casting  the  least  squares  problem  into  the  form  of  a  system  of  linear 
equations  and  then  utilizing  the  generalized  inverse  Pt  =  (P* P)-1  P*  to  solve  for  the 
coefficients. 


6.1  Example:  WS  Decomposition 

Example  1:  Heat  Equation  with  Dirichlet  Boundary  Control 

Consider  a  system  defined  by  the  partial  differential  equation, 

ft  =  6  ;  -’(0)  =  o,  -'(i)  =  «(<) 

l  y(t)  =  z\x0,t) 

This  system  has  transfer  function  [4] 

sinh  y/sx o 
Gi(s)  =  imhyr 

By  assuming  lowpass  characteristics  of  the  sensor,  we  can  write  the  overall  transfer 
function  as  , 

1  sinh  \/sxq 

Gr(s)  — - : — - — -p=- 

s  +  7r  sinh  y/s 

A  decomposition  of  G(s)  was  computed  using  12  dilation  levels  (m  =  —5, . . . ,  6)  and  up 
to  33  translations  at  each  dilation  level.  The  results  of  this  decomposition  are  shown  in 
Figures  6-8  which  are  different  representations  of  the  magnitude  of  the  wavelet  system 
expansion  coefficients.  Along  the  dilation  axis,  zero  corresponds  to  the  lowest  dilation 
level  (m  =  —5),  and  along  the  translation  axis  zero  corresponds  to  n  =  0.  Figures  6  and 


(41) 
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7  are  3D  and  contour  plots  respectively  of  the  coefficient  magnitudes  and  Figure  8  is  a 
density  plot  in  which  each  rectangle  is  shaded  according  to  the  corresponding  coefficient 
magnitude.  The  key  feature  in  this  decomposition,  which  is  probably  most  obvious 


Figure  6:  Wavelet  system  decomposition  of  heat  equation  transfer  function  -  3D  plot  of 
magnitude  of  expansion  coefficients 


Figure  7:  Wavelet  system  decomposition  of  heat  equation  transfer  function  -  contour 
plot  of  magnitude  of  expansion  coefficients 

in  the  density  plot  (Figure  8)  is  that  the  magnitudes  of  the  coefficients  are  very  well 
concentrated.  This  feature,  which  is  due  to  the  time-frequency  localization  properties  of 
affine  wavelets,  permits  us  to  pick  a  finite  number  of  ‘significant’  terms  in  the  wavelet 
system  decomposition.  In  doing  so,  a  finite-dimensional  approximation  to  the  original 
transfer  function  is  obtained.  Using  the  results  of  Section  5.3,  it  is  also  possible  to 
immediately  write  down  a  minimal  state  space  realization  of  the  approximating  system. 
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Figure  8:  Wavelet  system  decomposition  of  heat  equation  transfer  function  -  density 
plot  of  magnitude  of  expansion  coefficients 

7  Rational  WS  Approximations  to  Nonrational  H2(n+) 
Transfer  Functions 

Wavelet  system  decompositions  as  in  (16)  provide  a  means  of  representing  nonrational 
transfer  functions  in  H2(II+)  as  infinite  sums  of  rational  transfer  functions1.  As  men¬ 
tioned  earlier,  our  primary  objective  in  deriving  such  a  decomposition  is  to  devise  a 
systematic  method  of  constructing  rational  approximations.  What  is  required  now  is 
a  mechanism  which  allows  for  judicious  selection  of  a  finite  number  of  terms  from  the 
expansion  in  (16).  For  this  purpose  we  utilize  the  localization  properties  afforded  us  by 
affine  wavelets. 

For  a  significant  class  of  transfer  functions  arising  from  physical  systems,  the  WS 
decomposition  will  result  in  a  reasonably  compact  representation  in  the  sense  of  well 
localized  and  rapidly  decaying  coefficients.  This  phenomenon  can  be  explained  on  the 
basis  of  the  time-frequency  localization  properties  of  the  WS  transfer  functions  and  the 
observation  that  H2(fl+)  transfer  functions  arising  from  physical  systems  are  often  well 
localized  in  time-frequency  as  well.  It  is  possible  to  devise  a  variety  of  schemes,  each 
relying  on  time-frequency  localization  properties,  for  the  selection  of  terms  in  a  WS 
decomposition  for  use  in  a  finite-dimensional  approximation. 

One  possible  choice  for  the  selection  of  a  finite-dimensional  WS  approximation  can 
be  based  upon  knowledge  of  the  time-frequency  concentration  of  the  (nonrational)  trans¬ 
fer  function  which  is  to  be  approximated.  For  example  given  that  the  time-frequency 
concentration  of  a  transfer  function  /  €  H2(II+)  is  Q(f)  we  can  select  a  subset  of  the 
wavelet  system  transfer  functions  based  upon  the  size  of  Q(f)  f)  Qrn  n,  where  Qm  n  is  the 
concentration  of  the  WS  transfer  function  Gm’n(s).  Daubechies  in  [6]  provides  bounds 
for  the  error  of  such  approximations  in  terms  of  the  energy  of  /  outside  Q(f).  A  second 

1Note  that  causality  is  preserved  due  to  the  property  described  by  Equation  (9) 
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selection  scheme  could  consist  of  first  computing  the  wavelet  expansion  coefficients  for  a 
large  number  of  terms  and  then  simply  discarding  those  terms  with  ‘small’  coefficients. 
This  second  method  is  similar  to  those  used  to  achieve  signal  or  image  compression  via 
wavelet  decompositions  (c.f.  [3,  18]).  One  such  selection  criterion,  can  be  based  upon 
the  l 2  norm  of  the  coefficients.  Let,  us  assume  that  we  have  computed  the  WS  decom¬ 
position  of  a  transfer  function  G  and  let  am  n  =  (G,  S~lxV m<n)  denote  the  coefficients. 
We  choose  the  K  largest  coefficients,  (whose  indices  we  store  in  an  index  set  J),  where 
K  =  ff(ff)  is  the  smallest  integer  such  that, 


Ei 


E( 


i,n) 


a 


m,n  j 


where  0  <  8  <  1  is  some  predetermined  tolerance.  We  now  define  a  rational  approxi¬ 
mation,  Gj ,  to  G  as,  Gj(s)  =  E(m  n)eJ  Gm,n{s).  Note  that  the  results  of  Section  5.3 
immediately  provide  a  minimal  state  space  realization  of  Gj. 


7.1  Approximation  Error  Bounds 

The  following  lemma  provides  a  bound  on  the  approximation  error  using  the  scheme 
just  described. 

Lemma  7.1  Let  {xn}n  eZ  be  a  frame  for  a  Hilbert  space  Ft,  with  frame  bounds  A  and 
B.  Assume  ||a;n||  =  1,  for  all  n  £  Z.  For  any  f  £  ‘H,  define  an  approximation  f  to  f 

by, 

7=  ^2  if,  S~'xn) 

n£j 

where  J  is  an  index  set,  chosen  to  satisfy 

2l(/,s->x„)|2>(i-«)^!(/,s-'x.l)|2, 

neJ  neZ 


for  0  <  S  <  1.  Then, 

Proof:  By  the  frame  condition, 

b-,ii/-/ii2<E|(/-/"s‘1i*)|2 

JfceZ 


<  a-1]]/  -  7ll2- 


Since,  /  —  /  =  ffniJ  (/)  S~lxn)  xn,  we  Lave  two  coefficients  sequences  representing  the 
expansion  of  /  —  /  with  respect  to  the  frame  {xn}.  Therefore  by  Theorem  2.1, 

niJ  keZ 
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Therefore, 


fi-'ll/-/l!2  <  W(/-/,s-nt)| 

icX 

<  Y.[{f,S-'xn)\2 

n£J 

keT. 

<  M-1  ll/ll2 

from  which  the  result  follows.  ■ 

Remark:  The  error  bound  in  Lemma  7.1  is  established  in  the  general  setting  of  frames 
in  Hilbert  spaces.  When  applied  to  the  specific  case  of  affine  wavelet  frames,  the  bound 
may  prove  to  be  quite  conservative.  This  is  because  time-frequency  localization  proper¬ 
ties  of  affine  wavelets  are  not  exploited  in  the  lemma. 

7.2  Example:  Rational  WS  Approximation 

As  an  example,  consider  the  WS  decomposition  of  the  heat  equation  transfer  function 
in  Section  6.1.  Letting  8  —  0.4  the  above  described  scheme  results  in  the  selection 
of  7  terms,  with  a  corresponding  normalized  L2  approximation  error  of  0.109,  and  an 
approximating  system  of  dimension  22. 

8  Conclusions  and  Discussion 

In  this  paper  we  have  introduced  a  new  decomposition  (which  we  call  a  wavelet  system 
decomposition  )  of  H2(n+)  .  WS  decompositions  provide  a  representation  of  H2(II+) 
as  an  infinite  sum  of  rational  functions  of  bounded  degree.  Construction  of  the  WS 
representation  is  based  on  an  appropriate  grouping  of  terms  in  an  affine  wavelet  frame 
decomposition  of  H2(n+)  where  the  analyzing  wavelet  is  rational.  For  a  rational  an¬ 
alyzing  wavelet  of  degree  N  the  terms  in  a  WS  expansion  are  either  of  degree  N  or 

2  N. 

t 

WS  decomposition  naturally  leads  to  methods  for  constructing  rational  approxima¬ 
tions  to  nonrational  transfer  functions  in  H2(n+)  .  Criteria  for  selecting  a  finite  number 
of  terms  from  the  infinite  wavelet  system  expansion  rely  on  compactness  of  representa¬ 
tion  which  arises  due  to  time-frequency  localization  properties  of  the  wavelet  systems. 
Much  work  remains  to  be  done  to  understand  the  approach  of  the  present  paper  against 
the  background  of  prior  work  on  L2  and  H°°  approximation  theory  (c.f.  [11,  10,  12,  1]). 

A  point  should  be  made  regarding  our  motivation  for  constructing  affine  frames  in 
the  frequency  (Laplace)  domain  as  opposed  to  directly  decomposing  weighting  patterns 
via  affine  wavelets  in  the  time  domain.  First  of  all,  our  approach  naturally  preserves 
causality  of  the  approximating  system  since  each  term  in  the  WS  expansion  is  in  H2(n+) 
and  therefore  corresponds  to  the  transfer  function  of  a  causal  system.  If  one  were  to 
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use  affine  wavelets  in  the  time  domain,  special  ‘tricks’  need  to  be  applied  to  preserve 
causality  since  (even  in  the  case  of  compactly  supported  wavelets),  translations  would 
eventually  result  in  a  noncausal  wavelet.  Secondly,  there  would  be  no  mechanism  for  re¬ 
taining  rationality  of  the  Laplace  transforms  of  the  individual  wavelets,  since  translations 
(delays)  prevent  this. 

As  a  concluding  remark,  we  should  mention  that  the  following  result  has  recently 
come  to  our  attention  [16] 


Theorem  8.1  There  exist  no  orthonormal  bases  of  wavelets  for  H2(fl+)  generated  via 
the  framework  of  multiresolution  analyses,  with  an  analyzing  wavelet  xf  such  that  if  is 
continuous  and, 


<  c 


for  any  a  >  1/2. 


In  particular  this  rules  out  smooth  and  well  localized  orthonormal  bases  of  rational 
wavelets  for  H2(II+)  .  Although  it  has  not  been  demonstrated  that  multiresolution  anal¬ 
yses  are  the  only  means  of  constructing  wavelet  orthonormal  bases,  all  known  examples 
of  orthonormal  wavelet  bases  can  be  associated  with  multiresolution  analyses.  It  re¬ 
mains  an  open  problem  to  construct  an  orthonormal  wavelet  basis  (even  for  L2(IR)  ) 
which  does  not  arise  in  this  way.  The  above  theorem  suggests  that  if  one  wishes  to  use 
rational  analyzing  wavelets,  it  is  perhaps  necessary  to  consider  the  more  general  setting 
of  frames  (including  Riesz  bases). 
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Appendix  A 

Proof  of  Proposition  4.2 

Without  loss  of  generality  we  can  consider  the  case  m  =  0.  Therefore  we  need  to 
check  that  for  b  G  IR, 


/R 


f{u)g(uj  +  b)du 


-(/. 


f(u)g(u  -  b)du  )  =  0. 


(A.43) 


Consider  the  real  and  imaginary  parts  of  the  l.h.s.  of  Equation  (A.43)  separately. 
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Real  part  Taking  the  real  part  of  Equation  (A. 43)  , 


I  f{u)g(u  -  b)dJjj 


=  3?e 
=  Ue 


\j  ~ 

$Re  (J  f(u>)g(u>  +  b)du  -  J  f(u)g(u  -  b)du?j 

o  r°°  \ 

f(u)  (g(u  +  b)  —  g(u  —  b))  du  +  /  f(u)  (g(u  +  b)  —  g(u  —  b))  du  J 
-oo  J  0  / 

•0  r  oo  \ 

f(—u)(g(—u  +  b)—g(—u  —  b))du  +  J  f(u)(g(u  +  b)  —  g(u  —  b))duj 


=  Ke 


3?e 


/ 


f(—u)(g(—(u  —  b))  —  g(—(u  +  b)))du  +  j  f(u)(g(u  +  b)~g(u  —  b))du 

\ 

/(w)  (<7(0;  —  6)  —  </(u;  +  b))  du  +  /  f(u)(g(u  +  b)—g{u  —  b))du 


\  «/  U  ^ 

V  h(«) 

r  oo  \ 

[tf  (w)  -  T/fw)]  du  =  0 


Imaginary  part 

1m  f(<u)g(u  +  b)du  -  (^J  f(u)g(u-b)du^j 

=  Im  ^  J  f{co)g(u  +  b)du  +  J  f{u)g(u  -  b)du^j 

—  Zm  ^ J  f(u)(g(u-\-b)  +  g(u  —  b))du  +  J  f{u)(g(u  +  b)+g(u  —  b))diuj 
=  I m  J  f(—u)(g(—u  +  b)  +  g(—u  —  b))du  +  J  f(u>)(g(u  +  b)  +  g(io  —  b))dujj 

aoo  f  °° 

f(-u)(g(—(u-b))+g(-(u  +  b)))du  +  /  f{u)(g(u  +  b)+g(u 

(s 

yJ  o 


=  Im 


f(u)  ( g(u>  —  6)  +  g(u;  +  6))  (Zu;  +  J  f{u)  (g(u  +  b)  +  g(u  —  b))du 


%) 


Im  /  [tf(w)  +  tf  (w)]  dw  )  =  0 


Which  proves  the  proposition  ■. 


Appendix  B 

Proof  of  Theorem  5.1 
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We  outline  here  the  proof  of  Theorem  5.1. 


Proof  of  part  (a):  Since  (31-33)  is  a  2iV-dimensional  realization  of  Gm’n(s),  we  need 
only  prove  that  the  rational  function  Gm,71(s)  is  of  order  2N  for  all  n  6  2.  it  /  0  if  and 
only  if  (39)  holds. 

The  poles  of  Gm'n(s)  (see  (21))  are, 

n) ,  n) ,  vk(m ,  n) ,  I%(m,  n) } 

where  T]k(m,n )  and  i/k(m,n)  are  defined  by  (20).  Assuming  (39)  holds,  7]fc(m,n)  and 
vk(m,n)  have  nonzero  imaginary  parts  for  all  m,n  £  Z.  Therefore,  at  s  =  r]k(m,n)  for 
example,  the  term, 

~  Vk{m,  n))(5  _  n)) 

k 

which  appears  in  the  numerator  polynomial  Nmn[s )  of  Gm<n  will  be  nonzero.  Further¬ 
more,  the  assumption  that  T  is  of  order  N,  implies  that  pj  ^  zk  for  all  z,  k.  Therefore 
at  s  =  the  remaining  (nonzero)  term  in  the  numerator  of  Gm'n  will  be, 

Nm>n(7?fc(m,ra))  =  a  JJ(s  -  /9J(m,n))(s  -  JJ(s  -  rj]:(m,n))(s  -  Vk(m,n)), 

j  k 

which  shows  that  Gm'n  has  no  zeros  at  rjk(m,  ri).  Exactly  the  same  argument  can  be 
used  to  show  that  there  is  no  cancellation  with  zeros  for  the  remaining  poles.  Thus  Gm’n 
is  of  order  2 N. 

The  converse  is  readily  proved  using  a  similar  pole-zero  cancellation  argument.  H 


Proof  of  part  (b):  By  hypothesis,  {Am  n,  Bm  n,  Cm<„)  is  a  minimal  realization  of 
Gm’n(s)  =  Nm>n(s) / Dmtn(s)  for  all  (m,  n)  E  J.  Therefore  Nm>n  and  Dm  n  are  coprime 
polynomials.  Consider  the  parallel  combination  of  two  WS  transfer  functions. 


Gj{s)  =  G 

m\  ,7ii  (s)  +  G 

m2  ,n2  (*), 


(B.44) 


where  J  ={(mi,  nj),  (ra2,  n2)}.  In  this  case 

{Nmi,ni(s)i  Dm2,n2{s))  and  {Nm2tn2(s)Dmuni(s)),  are  relatively  coprime  pairs.  Rewrit¬ 
ing  (B.44), 


Gj(s)  = 


Nj(s) 

Dj(s) 

Nmiini{.s)Dm2tn2(s)  +  Nm2^(s)Dmitni(s) 
Dm\,n\  iS)^m2,n2  (5) 


Therefore  Nj  and  Dj(s)  will  be  relatively  coprime  if  and  only  if  Dmi  Ul  and  Dm2<n2(s) 
are  relatively  coprime.  The  roots  of  Dmuni  and  Dm2  ,„,(s)  are  given  by  (21).  Thus 
and  Dm2  n2(s)  will  have  a  common  root  if  and  only  if, 


ao  mi  (Pi  +  =  «o  m2{Pj  ±  innb0), 
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for  some  /  andj.  Solving  for  nl, 

°o 

=  — —  (a™1_?ri2(3fe  pj  +  Tim  pf)  —  ('lie  pl  +  iXm  pfj)  ±  a™l~m2??.2 

—  (a™1-™2  I’m  pj  -  1m  pt)  ±  «™1_ni2ra2 

4Kmrra,&pr3fep;)  (B.45) 

However,  by  hypothesis  (40)  there  can  be  no  integer  solution  to  (B.45)  since  the 
imaginary  part  of  (B.45)  will  always  be  nonzero.  The  above  argument  for  two  WS 
transfer  functions  is  easily  extended  to  any  finite  number  of  WS  transfer  functions  in 
parallel.  ■ 


Appendix  C 

Computation  of  Frame  Parameters 

Given  an  admissible  analyzing  wavelet  g ,  the  following  theorem  of  Daubechies  is 
useful  in  determining  dilation  step-size  a  and  translation  step-size  b  for  which  (</,  a,  b) 
generates  an  affine  frame  for  L2(1R)  . 

Theorem  C.l  (Daubechies  [6])  Let  g  €  L2(1R)  and  a  >  1  be  such  that: 


(i) 

m{g ;  a)  =  ess  inf  Y  \g(anu))\2  >  0 

|u/|€[l,a]  1 — ' 
n 

(C.46) 

(2) 

M{g\a)  =  ess  sup  5  \g(anu)\2  <  oo 

(C.47) 

Me[i,«]  „  “ 

(S) 

OO 

lim2  /?(27rfc/6)1/2/5(  —  2irk/by/2  =  0, 

(C.48) 

ffis)  =  ess  sup  Y'  \g{anu)\\g(anu>  -  s)|. 

Then  there  exists  Bc  >  0  such  that  (g,a,  b)  generates  an  affine  frame  for  each  0  <  b  <  Bc. 


Proof  of  the  following  corollary,  can  also  be  found  in  [6]. 
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Corollary  C.l  If  g  £  L2(IR)  and  a  >  1  satisfy  the  hypotheses  of  Theorem  C.l  then, 


Bc  >  bc  =  inf  { fo|  m(g-  a)  —  2  ^  (3(27rk/by/2[3(—2irk/by/2  <  0} 

k= 1 

mid  for  0  <  b  <  bc,  the  frame  bounds  A  and  B  can  be  estimated  as, 

OO 

A  >  b-\m{g-a)-2S£jP{2irklbyl2(d{-2^klbyl2) 

k=l 

OO 

B  <  b-1{M{g-,a)  +  2j2^i^k/byf-2y(-2Trk/b)1B) 

k~l 


(C.49) 


(C.50) 


26 


